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SUMMARY 


A method is presented for determining the lateral -stability deriva- 
tives, tr ans fer-function coefficients, and the modes for lateral motion 
from frequency-response data for a rigid aircraft. The method is based 
on the application of the vector technique to the equations of lateral 
motion, so that the three equations of lateral motion can be separated 
into six equations. The method of least squares is then applied to the 
data for each of these equations to yield the coefficients of the equa- 
tions of lateral motion from which the lateral-stability derivatives and 
lateral -motion transfer-function coefficients are computed. Two numeri- 
cal examples are given to demonstrate the use of the method. 


INTRODUCTION 


In the reduction and generalization of flight-test data, whether for 
loads, stability, or control purposes, the airplane stability derivatives 
and the coefficients of the transfer functions are often required. A 
great deal of emphasis, therefore, has been placed on the development of 
analytical methods for reducing flight data to obtain these basic deriva- 
tives and coefficients. 

A number of recent methods, for example references 1 to b, are now 
available for analyzing longitudinal maneuvers and determining the 
longitudinal-stability derivatives and transfer-function coefficients 
from flight data. References 1 and 2 present methods of determining the 
longitudinal -stability derivatives and transfer functions directly from 
transient data. Reference 3 reduces data for longitudinal motion deter- 
mined from the forced-oscillation technique by means of circle diagrams 
to longitudinal-stability derivatives and frequency response. Mueller, 
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in reference 5 ) was one of tlie earliest to use vector representation in 
the equations of longitudinal motion to represent the derivatives and 
integrals of the variables. Schumacher, reference 4, represented the 
frequency responses to longitudinal motion as vectors and substituted 
them into the equations of longitudinal motion and the transfer func- 
tions. He then applied the method of least squares to these vector 
equations, a method which he found very effective in deter mining certain 
of the longitudinal -stability derivatives and. transfer-function 
coefficients . 

The problem of analyzing lateral motions, however, has not received 
the same amount of attention as that for longitudinal motion, perhaps 
because it is more complicated. Several analytical investigations have 
been undertaken and a few methods have been proposed such as the circle- 
diagram method (ref. 6 ), the step-function-response method (ref. 6 ) in 
which the response of an airplane to a step deflection of the rudder or 
aileron is analyzed, and the free -oscillation method (ref. 7 ) in which 
the period and damping of the free vibrations of the aircraft due to a 
pulse-type input are analyzed. Since the usefulness of these methods 
is limited by the number of derivatives which can be extracted, there is 
still a need for a more general method of analy sis that will extract all 
the significant lateral-stability derivatives from flight data. 

It is the purpose of this paper to present a method for determining 
the lateral-stability derivatives of a rigid airplane and to illustrate 
its use by applying it to two examples. The method is based on the 
vector representation of the frequency responses to lateral motions. 

This vector approach permits separation of each of the equations of 
lateral motion into a real and imaginary equation. A least-squares 
method is then applied to the data in each of these equations or combi- 
nations thereof to yield the coefficients of the equations of lateral 
motion. The lateral-stability derivatives and transfer -function coef- 
ficients are then determined from these coefficients and the known air- 
craft mass parameters. 

The method is applied to two specific examples, one in which the 
frequency responses to a rudder input are known and one in which the 
transient responses to aileron deflection are known. In the latter case, 
the frequency responses were obtained from the transient motions by two 
methods and the stability derivatives computed. 

An attempt has been made to schedule the procedure so as to reduce 
the dependence of the results obtained from this method on the deriva- 
tives that can be least accurately obtained from the particular data 
being analyzed^ however, further improvements may be made as further 
experience is gained in the application of the method. 
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SYMBOLS 


®y 


Ajf^ 

B P> B cp> B ay 


b 

C L 



lateral acceleration, ft/ sec 2 
parameters defined by equations (8) 

parameters defined by equations ( 9 ) 

wing span, ft 

lift coefficient, L/qS. 


Cn 

c* 

c Y 

Cl, c 2 , c^. 


D 


F l> F 2' F 3 
G 2 , Gj 


g 


Iz 


yawing-moment coefficient, il/qSb 

rolling-moment coefficient, L'/qSb 

lateral -force coefficient. Lateral force/qS 

coefficient of transfer functions 
(defined in table H) 

differential operator, ~ 

dt 

natural logarithmic base 

forcing -function coefficients representing rudder 
effectiveness (defined in table I) 

forcing-function coefficients representing aileron 
effectiveness (defined in table I) 

acceleration due to gravity, ft/sec 2 

moment of inertia about stability X-axis, inky 2 , 
slug-ft 2 

moment of inertia about stability Z-axis, mkg 2 , 
slug-ft 2 
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I XZ 

i 

%, k 2 , k 3 , 

k x 

% 

KZ 

K xz 

k Xo 

L 

L* 

M 

m 

N 

p or 0 


product of inertia referred to stability axes (nega- 
tive when the positive direction of the X principal 
axis is inclined above the flight path, i.e., when 
q is positive) 

stability-derivative coefficients of the equations 
of lateral motion (defined in table I) 

radius of gyration about stability X-axis, ft 

radius of gyration about stability Z-axis, ft 

nondimens ional radius of gyration about longitudinal 
stability axis, 

nondimens ional radius of gyration about vertical 
stability axis, 

nondimens ional product-of -inertia parameter, 

- IK/*) 

radius of gyration about principal longitudinal 
axis, ft 

radius of gyration about principal vertical 
axis, ft 

lift, lb 

rolling moment, ft -lb 

Mach number 

aircraft mass, W/g, slugs 


- (%/>) 2 


cos q sin q 





yawing moment, ft -lb 

incremental roll in g angular velocity about 
X-axis, radians/sec 
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q. 

r or t 

R 

S 

s 

t 

V 

V 
W 

X, Y, Z 
a 

P 

5 a 

8 r 

n 


Pb 

p 

T 

$ 

0 

6Q 


dynamic pressure, ipV*^, fb/sq. ft 

incremental yawing angular velocity about 
Z-axis, radians/sec 

amplitude ratio 

wing area, sq. ft 

Laplace transform variable 

time, sec 

incremental component of velocity V along 
the Y-axis, ft/sec 

true airspeed, ft/sec 

aircraft weight, lb 

airplane stability axes (see fig. 1) 

angle of attack, radians 

angle of sideslip, v/V, radians 

aileron control deflection, radians 

rudder control deflection, radians 

inclination of principal longitudinal axis of inertia 
with respect to flight path (positive when the posi- 
tive direction of the X principal axis is inclined 
above the flight path) 

nondimens ional mass parameter used for lateral 
equations, m/pSb 

air density, slugs/cu ft 

time parameter, m/pSV, sec 

phase angle, radians 

angle of roll, radians 

angle of yaw, radians 

angular frequency, radians/sec 
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Bar notation: 

a bar over letter represents ma-xi mum value 

j a J bars on sides of symbol represent absolute value 

Matrix notation: 


[] 

{} 


square matrix 


column matrix 


The lateral-stability derivatives are expressed by subscript nota- 


tion as, for example: = 




3 83 


Jn P 


— *^n 


<£> 

2V 


and. C 


“r 


^n 

W 

2V 


Phase angles are also indicated by subscript notation as <J> D is 

P 5r 

the phase angle between the input rudder deflection and the output side- 
slip angle. 


DEVELOPMENT OF METHOD 
Equations of Lateral Motion 


The equations of lateral motion based on the stability axes (as shown 
in fig. l) and on the usual assumptions of linearity, smal 1 angles, and 
maneuvers which start from a level-flight condition are usually written 
in the following form: 

(a) For the sideslipping motion, 

(D + K-l) 3 - + Dt = F-^Ct) (r) 

or, since ay = v(d£ + Di|r - K 2 0), 

iay + Kq3 = FpS^t) (la) 

(b) For the rolling motion, 

Kj3 + (d 2 + K^d) 0 - (f^D 2 + KgD) * = F^t) + G^Jt) 


( 2 ) 
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(c) For the yawing motion, 

-Y - + v) ^ + + he?) * = F 3 b r M + G 3 5 a (t) (5) 


where the K, F, and G coefficients are defined in table I. The 
forcing functions of these equations are written for rudder inputs and 
aileron inputs combined; for the sake of simplicity in developing and 
demonstrating the method, only rudder inputs F^6 r (t), Fgb^t), 

and FjS r (t) are considered. The method, however, applies equally well 

to both aileron and rudder inputs as is shown subsequently. 


Vector Interpretation of Frequency Response 

In many current studies of airplane dynamics, the frequency response 
of the aircraft is determined from flight data. Having the data avail- 
able in the frequency plane offers certain advantages over having the 
data in the time plane, since use can then be made of the vector inter- 
pretation of frequency response. The development which follows will 
assume that the data to be analyzed are available in frequency-response 
form. In the event that only transient-response data are available, the 
transformation to the frequency plane can be made from a selection of one 
of the several methods compared in reference 8. 

The vector technique is applied to the determination of lateral- 
stability derivatives in the following manner . If an input ^ to a 
linear system (a system described by a linear differential equation) is 
considered to have a sinusoidal variation of frequency cd, it can be con- 
sidered to be a vector of -unit magnitude lying along the real axis of the 
complex plane as shown in figure 2 (a). The first derivative of the sinus- 
oidal input D6r is then obtained by multiplying the amplitude by 00 and 
rotating the resulting vector 90° counterclockwise in the complex plane 
which is equivalent to advancing the phase angle by 90 °. Each succeeding 
higher derivative is found by multiplying by to and rotating the vector 
90° in the complex plane; thus. 


H 

U 

w 

DBj. = io 

( 5 ) 

D^S-r = -afi 

(6) 
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The steady-state output or response of the linear system to a sinus- 
oidal input may he considered to he a vector of magnitude equal to an 
amplitude ratio R and having a direction angularly displaced from the 
real axis hy a phase angle $>. Thus, in the complex plane, the vector 
for sideslip may he represented as 


3 = 


JL 

Sr 


i$c 


Sr _ 


JL 

Sr 


^cos + i sin = Ap + iBp 


By definition. 


Ap = 
— 


B 3 S 




Sr 

± 

§r 


Sr 

hr 


cos <5 1 


Sr 


cos 


3>\)r 


&r 


sin <£ 


3* 


sin <5 




Sr 






£ 


COS 


Sr 


■^y 


B 0- 


Sr 


cos G> 




P 




sin 


05. 


B„ = 

®y 


fy 

hr 


sin $ 


a y, 


hr 


The derivatives of p are represented as 


( 7 ) 


( 8 ) 


( 9 ) 


Dp = CD 


= CD 


JL 

Sr 

JL 

Sr 




-sin + i cos <&p^ 


= -dBn + 


p + iaA p 


( 10 ) 
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and 

D 2 p = ^-Le i f P ^ + *) 

Sr 


- ^ 1 ; (-° os V - 1 sln V) 


= -trpAp - i<r?Bp 

(H> 

Similarly, for roll (see fig. 2(h)), the equations corresponding to equa- 
tions ( 7 ), ( 10 ), and (ll) would he 

$ = ^ + “0 

( 12 ) 

D$ = -aiB^ + l(xA$ 

(13) 

= kb 2 A^ - io> 2 B 0 

(14) 

for yaw. 


t = Ay + iBy 

(15) 

Dijr = -city + icoA^ 

( 16 ) 

D 2 ^ = - io^B^ 

(17) 

and, for lateral acceleration. 


ay = Aa^ + iB^ 

( 18 ) 


Substituting equations (4) to ( 18 ) into equations (l), (la), ( 2 ), and ( 3 ) 
»nri equat ing the real and imaginary values in each equation yields the 
following eight simplified equations: 
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From equation (l) : 


Kj_Bp - ( a ^ + A,^ (19) 

and 

K x Ap - = F x + + B^) (20) 

From equation (la) : 

VKiBp = -B^ (21) 

VK-jAp = VF X - (22) 

From equation (2) : 

K^Bp + K],aA0 + K^B^ - KgcoA^ = o£b0 (23) 

and 

KjAp - + K^A^ + - F 2 = a?Af (2k) 

From equation (3): 


-KyBp + Kqo£b 0 - K^oA^ + K 10 oA^ = 

and 

-Ky-Ap + Kqcj3 2 A^ + K^oiB^ - K 10 aA|; - *3 = 03 2 A^ 


(25) 


(26) 


In equations (19) to (2 6) the 


8114 B p,^,+,ay 


terms will 


he available at the particular values of cd from the frequency-response 
curves. It will he assumed that the Kq term is determined from the 
velocity measurement and the and Kg terms, which are equal to the 

ratios lyz/ly ^ I-^/l^, respectively, are known from either weight 

and balance calculations or measurements so that only the values 
of Kq_ 3 4 g y ^ io H - r1 ^ ^1,2,3 8X6 "^° Be determined from the data. 

(it should he noted that even though K = and Kg are assumed to he known 


a need still exists for a simple method of determining aircraft moments 
of inertia and the location of the principal longitudinal axis.) 
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Equations (19) to (2 6) are now used to compute the K and F coef- 
ficients. The advantage of separating each of the equations of lateral 
motion into real and imaginary equations now becomes apparent. The two 
derivatives most difficult to determine are C 7 and C„ cont aine d in 

i r Up 

the Kg and Ky coefficients, respectively. The Kg coefficient may 
be eliminated between equations (23) and ( 2 k ); similarly the K^ coef- 
ficient may be eliminated between equations (25) and (26). The equations 
thus obtained, when used to compute the unknown coefficients, result in 
better conditioned matrices and more accurate results. In the us ual case 
it has been found u nn ecessary to perform the elimination of the 
K^ coefficient . 

As a first step in solving for the K and F coefficients, equa- 
tions (19) to (26) are fitted to the frequency-response data over a range 
of values of 0 by the method of least squares . The theory of least 
squares is well-known and is derived in many textbooks (e.g., ref. 9) and 
no attempt is made here to repeat its derivation. 

The application of the least-squares method to equations (19) to (26) 
converts these equations to their computational form. It is recommended 
that, if the lateral-acceleration frequency response is available, equa- 
tions (21) and (22) be used in preference to equations (19) and (20) for 
the computation of Kq and Fq. 

Equation (19) may be expressed simply as 



(27) 


and yields Kq. Similarly, equation (20) can be expressed as 


F 1 = 5 (kqAp - - niBp - 


(28) 


and yields Fq. 
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For the same variables K^_ and *L, the least -squares forms of 
equations (21) and (22) for use with the lateral-acceleration frequency 
response are: 



Equation (25) can be expressed as 





P -'X 

Kt 


£ (A 


l^[H3 2 



> «= < 


|WiNi 


^[^2 2 _ 

K 10 




and yields K^., K^, and K-^q. 

Equation (26) can be expressed as 


*3 = i 5_ (-Vp + * 8 “% + W ‘ K io^ - m \) - (52) 

and yields . 

In order to evaluate the remaining constants, K^, Kg, and f 2 , 

the coefficient Kg is first eliminated between equations (23) and (2k). 
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This step is taken in order to obtain a more accurate value of the coef- 
ficient Kg, wMch. is less significant under these conditions, and to 
provide a better conditioned matrix of coefficients in the computation 
of K 4 , and Fg. Upon elimination of the Kg term between equa- 

tions (23) and (2k), the following relation is obtained: 


K^BpIty + ApA^ + - B^A^ - FgA^ = + A^A^ - 


± 
«r| 
(33) 


Upon the application of the method of least squares, equation (33) 
becomes 



|^(“¥V Vf)j ^ l (-A)j( B P I V + 

% 

ft (»& + t 

J-fflL 


- K)p 



F 2 



which yields K^, K^, and Fg. 

The Kg coefficient may now be obtained from equation (24) which, 
after the method of least squares has been applied, may be expressed as : 


*6 = 



+ • oP A^ - Kjp<iP 


Ki^oiB^ - KjA^ (afi^ 



(35) 
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An alternate expression for Kg may be obtained from equation (23) 
which may be expressed as 


*6 = 


<5=<qL V 3 


£ IH «T 


(36) 


Values of Kg, which define the derivative C^, are Important in 

the spiral mode. In order to obtain accurate values of Kg, aircraft 

motions which bring out the full effect of the long-period spiral mode 
should be analyzed to obtain accurate frequency-response data near oo = 0. 


Tn the case where the solution for 



is found to involve ill- 


conditioned matrices, K^ should be eliminated between the real and 
imaginary equations (25) and (2 6). 


With all the K and F coefficients of the equations of motion 
determined, the lateral-stability derivatives can be computed from the 
K and F coefficients and the known aerodynamic parameters by use of 
the definitions of table I. The transfer-function coefficients can be 
calculated from the definitions in table H and the transfer functions 
and modes of lateral motion can be obtained from the equations in 
appendix A. 


SUGGESTED PROCEDURE 


In order to aid in the application of this method to the analysis 
of flight data, a suggested step-by-step procedure which is presented in 
this section has been worked out. An effort has been made to schedule 
the procedure so as to reduce the dependence of the results on the deriva- 
tives that can be least accurately obtained from the particular data being 
analyzed. Alternate steps are suggested where it was found that some par- 
ticular derivative might be more accurately determined by one or the other 
of two approaches under certain conditions. Wo weighting of specific 
groups of data is employed in the least-squares procedure j weighting can, 
however, be employed when it is considered desirable to put more dependence 
on data regarded as more reliable. 
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As a further demonstration of the method, two numerical examples 
have been carried out according to this procedure. In one example, pre- 
sented in appendix B, the lateral-stability derivatives are calculated for 
an airplane whose transient responses to aileron deflection are assumed 
to be known. This example was chosen in order to gain some insight into 
the dependence of the method on the accuracy with which the frequency 
response can be obtained from the transient response. 


The other example is for the hypothetical rigid airplane whose mass 
and geometric parameters are listed in table HI. The assumed frequency- 
response data of sideslip angle {3, roll angle f, yaw angle i|r, and 
lateral acceleration ay for a rudder-deflection input such as might be 
obtained from analysis of flight data are listed in table TV and plotted 
in figures 3 bo 6. The analysis of these data has been used to illu strate 
the following suggested procedure. 


(l) Tabulate parameters and working equations .- Tabulate the air- 
plane parameters and any stability derivatives that are known or can be 
computed directly from those that are known as illustrated in table HI. 
Tabulate the least-squares equations (27) to (32) and (3^) to (36) with 
any known terms on the right-hand side as has been done in equations (27a) 
to (32a) and (3^a) to (36a) in table V for the example. ^Since was 

taken to be zero in the example, all terms containing - were dropped; 
however, in the example of appendix B, was included and the 

term was determined.^ 


(2) Tabulate frequency-response data .- Tabulate the amplitudes and 
phase angles of the four lateral variables at the values of cd to be 
used in the analysis. For simplicity in the example, as shown in table IV, 
10 integral values of cd, evenly distributed over the range of the data, 
were chosen; however, in cases where more accurate points occur at non-, 
integral values of cd, it is advisable to use such points, where possible. 
The range of values of cd should be restricted to the rigid response. 


but the point cd 


0 


should be avoided because 



is infinite . 


(3) Compute vector components of data .- Compute the vector compo- 
nents from the data listed in step 2 and tabulate them as shown in 
table IV. It is advisable to use at least ^ --place tables of sine and 
cosine values in this computation. Additional columns of cdA^, oiBp, 

cdA^, . . . and crf^Ap, cii^Bp, a£/y, . . . will be useful in subsequent 
operations but are not listed for the example. 
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(4) Compute Kq.- Compute Kq from equation (27a). (See 

table Vl(a).) In this step and in subsequent steps, summation of tbe 
corresponding elements of tbe indicated products of tbe data columns 
evaluates tbe matrices of coefficients and knowns. Tbe resulting simul- 
taneous equations are readily solved by any standard method. 

(4a) Alternate step 4 .- Compute Kq from equation (29a). (See 
table Vl(c) . ) 

(3) Confute Fq .- Compute Fq from equation (28a) with tbe values 
of Kq from step (4). (See table Vl(b).) 

(3a) Alternate step 3 - - Compute Fq from equation (30a) with tbe 
values of Kq from step (4a). (See table Vl(d).) 

(6) Compute Ky and KqQ.- Compute Ky and Kqo from equation (31a) 
(See table Vl(e) . ) 

(7) Compute F^. - Compute F^ from equation (32a) witb tbe values 
of Ky and Kq 0 from step (6). (See table Vl(f).) 

(8) Compute K^, and F 2 .- Compute K^, K^, and F 2 from 

equation (34a) . Because of tbe combination of two equations to obtain 
equation (34) in tbe method, tbe combination of data columns is more 
complicated in this step than in others. (See table VI (g).) 

(9) Compute Kg.- Compute Kg by one of tbe following methods: 

(a) If F 2 ha s not been eliminated, calculate Kg from equa- 
tion (35a) witb the values of Kj, Kij., and F 2 from step (8). (See 
table VI (h) . ) 

(b) If f 2 has been e limina ted calculate Kg from equation (36a). 
(See table Vl(i).) 

(10) Compute stability derivatives .- Compute tbe lateral-stability 
derivatives from tbe relationships in table I and tbe now complete set of 
values of K and F. (The computed and known values of tbe K 

and F coefficients and tbe stability derivatives for tbe example are 
given in tables VH and VUE, respectively.) 
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(11) Compute transfer-function coefficients .- Compute the transfer- 
function coefficients from the relationships given in table II and the 
values of K and F. (The computed and known values of the transfer- 
function coefficients for the example are given in table IX.) 

(12) Compute the modes . - The modes of lateral motion may now be 
computed, if desired, by use of the equations in appendix A. (The values 
of the lateral modes for the example are given in table X. ) 

(13) Accuracy check .- When the true derivatives are unknown, the 
check on the computations is obtained in the following manner; Calculate 
the frequency response from the equations in appendix A and the computed 
coefficients of step (ll) . Compare these values with the data of 

step ( 2 ). 


DISCUSSION 


Data requirements .- The method of this paper is based on the assump- 
tion that equations (l) , ( 2 ), and ( 3 ) represent the lateral motions of the 
airplane. These linear differential equations are satisfied by small 
variations in roll, yaw, and sideslip but will not necessarily hold for 
large changes in these variables of in regions of flight where nonlineari- 
ties occur. It is therefore desirable to obtain data for transformation 
to the frequency plane from small motions recorded within ranges of Mach : 
number and other variables in which the coefficients are constant; also, ■ 
the range of frequency-response data to be used must be limited to values 
for which the aircraft acts as a rigid body and higher frequency instru- 
ment inaccuracies are avoided. Thus, the assumptions of the method will 
be most nearly in accord with the true physical nature of the problem. 

In any case, experience has shown that the derivatives obtained by this 
method of analysis will accurately reproduce the frequency response frem 
which they were extracted. 

Since the present method requires flight data in frequency-response 
form, either the forced-oscillation technique or the technique of oscil- 
lating the rudder or aileron at various frequencies and measuring the 
steady-state response in p, D f, D|r, and a y must be employed or the 

transformation of a transient response from the time plane to the fre- 
quency plane must be made by an appropriate method. A comprehensive sur- 
vey of methods for effecting this transformation from time to frequency 
plane is given in reference 8 which evaluates the Fourier integral method 
as well as the use of devices such as the Fourier synthesizer, Coradi 
harmonic analyzer, and IBM machines for obtaining frequency response from 
transient -response data. 
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If, as is usually the case, the transient responses in roll and yaw 
are measured in the form of time histories of D 0 and D\|r rather than 
of 0 and ijr, it is recommended that these time histories of D 0 
and Ity he converted to frequency response initially and then the fre- 
quency responses of (f> and he computed from the frequency responses 
of D0 and Ity hy use of the relations 


and 


D0 

= CD 

i 


Dir 

= CD 

JL 

§r 


§r 


§r 


§r 


x - x + 90 ° X “ X + 900 


Flight data are usually measured about the body axes; before they 
can be used, however, either they must be converted to a form corre- 
sponding to the stability axes or the development must be restated in 
terms of body axes. The following equations are presented to be used 
in transforming to the stability axes: 


0 = cos a + tt) sin a 

• • \ 

i|r = ^ cos a - sin a 


where a is the angle of attack (angle between the body X-axis and the 
stability X-axis) and the subscript b refers to the body axes. In 
addition, the sideslip-vane measurements can be converted to the corre- 
sponding values at the center of gravity by means of the formula 




cos a 


where the subscript v refers to the vane and is the distance from 

the center of gravity to the sideslip vane measured along the body X-axis. 

If the lateral acceleration By as recorded by an accelerometer is 

used, it must be corrected for any components of angular velocity and 
acceleration present in such records because of the displacement of the 
accelerometer from the airplane center of gravity. 
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Sensitivity of the method and possible simplifications . - Even though 
the foregoing general rules are observed, the question arises as to how 
the derivatives obtained compare with the actual derivatives of the air- 
plane, and how sensitive these derivatives are to inaccuracies in the x 
frequency-response data. These inaccuracies stem from several sources 
such as measurement of the transient motions, transformation from time to' 
frequency plane, and reading errors. Care should he exercised to keep 
these Inaccuracies to a m inim um. 


No attempt is made to present a comprehensive error analysis in this 
paper, hut several general remarks appear to be in order. A comparison 
of the values given in tables VII, VUE, and IX for the coefficients of 
the equations of lateral motion, the lateral-stability derivatives, and 
the transfer-function coefficients indicates essentially no errors due 
to the method itself. It has been found, however, that the derived values 
of Cy^ and Cy are sensitive to slight variations in the quanti- 
ties Ap, Bp, A^, and B^. (see eqs. (27) and (28)). This sensitivity 

i3 due to the fact that the airplane angle of yaw is very nearly equal 
to the negative of the sideslip angle and hence accurate determination 
of the lateral acceleration from measurements of ijr and (3 is diffi- 
cult. This difficulty can he obviated by a direct measurement of the 
lateral acceleration with an accelerometer and more accurate values 


determined for Cv and 


C Y 5 by means of equations (29) and (30), 


As brought out in appendix B, the accuracy with which the values of 
the derivatives C 7 . Cb, , and Cy,, can be determined from frequency- 

± J ’P 

response data appears to be dependent upon the inclusion of data for fre- 
quencies near the natural frequency of the airplane. This fact is brought 
out further by some unpublished results that indicate that derivatives 


such as C n and Cy 


are very sensitive to variations of the frequency 

In the case of Cv, the 


responses near the airplane natural frequency. ■LU UJJLW 

conclusion is obvious since this derivative is known to contribute largely 
to the airplane Dutch roll damping and it is also known that the peak value 


of the frequency-response amplitude ratio, say 


is highly dependent 


upon the airplane Dutch roll damping. It appears therefore that an 
accurate representation of the frequency responses at frequencies near 
the airplane natural frequency (for instance +50 percent) is probably 
very important in the determination of the lateral-stability derivatives. 


In general, when the number of unknown derivatives can be reduced, 
the remaining derivatives can be determined with greater reliability. 

For this reason, several possible simplifications are presented for the 
calculation of the derivatives C , C , C, , and C, . The airplane 

P S r P 5 a 
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natural frequency ±6, for Rmal 1 values of the product of inertia, well 
approximated "by the expression 

o qSb 

« Cnp — 

If the natural frequency ccy is apparent from the available transients 
and the inertia parameter is known, a very good value of can gen- 

erally he obtained frcrn this expression. 


Extrapolation of the amplitude ratio 


St- 


afford a means of obtaining the derivative C 


to zero frequency may 
subsequent to deter- 


n 6r 

mination of from the preceding frequency relation. The value of 


_P_ 

Sr 


at a> = 0 obtained by this extrapolation is not indicative of the 


actual static sensitivity, hut rather of the apparent static sensitivity 
which includes the very lightly damped spiral mode. This apparent steady- 

C nc 

state response is approximately equal to the ratio -31, from which 

% 


■'n. 


Sr 


can he estimated if C 




is known. 


The transfer function which relates the rolling velocity to aileron 
deflection is well approximated (particularly at low frequencies) by the 
expression 


i-do) 


G/- 


°2 


-i tan-1 r~ 


K 4 


Ki. + icD 


(K^ + cxf~J 


1/2 


Therefore, $$ (cd) = -tan"l — from which Kh , and subsequently Cj , 
» On K|^ 


5a 

can be determined. Also, by extrapolating the amplitude ratio 


L 


to 


cd = 0, the parameter Gg and hence can be obtained once Kj^ has 

^a 

I -LI 


been determined. As was the case for 


yield the true value of 


, this extrapolation does not 


at cd = 0, which is actually zero, but gives 


the apparent value due to the lightly damped spiral mode. 
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CONCLUDING REMARKS 


A method has been presented for extracting lateral-stability deriva- 
tives from frequency-response data which have been derived from aircraft 
transient responses to arbitrary control inputs. In order to demonstrate 
the use of the method, the lateral-stability derivatives have been calcu- 
lated for two hypothetical airplanes for which the frequency-response data 
or transient-response flight data were assumed to be known. Simplifica- 
tions were proposed for obtaining certain of the derivatives, and although 
no error analysis is presented in this paper, seme general observations 
can be made concerning the sensitivity of the method to inaccuracies in 
the original data. 

In view of the limited experience in the determination of lateral- 
stability derivatives from flight data, and particularly in using the 
present method, it appears that an investigation should be made to deter- 
mine the effects of some of the errors inherent in analyses of transient 
responses, such as recording accuracy, reading accuracy, and so forth, 
on the derived frequency responses and hence on the computed lateral- 
stability derivatives. It is noted here that, although the methods pre- 
sented in this paper for the extraction of lateral-stability derivatives 
from frequency-response data are theoretically correct and mathematically 
sound, the accuracy with which the lateral -stability derivatives of an 
airplane can be determined by this method is directly dependent upon the 
accuracy with which the frequency responses of the airplane are known. 
Also, some consideration should be given to the types of inputs likely 
to afford good frequency-response data and to testing procedures whereby 
some of the stability derivatives that are most difficult to determine 
may be accurately obtained. In addition, accurate frequency-response 
data for specific ranges of frequency may yield more accurately certain 
of the stability derivatives and this possibility should be investigated. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., February 1 6 , 195^ • 
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APPENDIX A 

DEVELOPMENT OF EQUATIONS FOR TRANSFER FUNCTIONS, MODES, 
AND FREQUENCY RESPONSE FOR LATERAL MOTION 


In this appendix the transfer functions are derived from the equa- 
tions of lateral motion and the modes and frequency response are expressed 
in terms of the tr ans fer-function coefficients . 


Equations of Motion 

The motions of a rigid aircraft resulting from a rudder-deflection 
input are assumed to he expressed hy the following three standard line- 
arized equations of lateral motion: 


(d + K-A p - + D<r = F 1 6 r (t) 

(Al) 

K^p + (d 2 + KjjD) <f) -(K^D 2 + K£D)* = F 2 5 r (t) 

(A 2) 

-I^p -(kqD 2 + I^d) 0 + (d 2 + K 10 d)+ = F 5 S r (t) 

(A3) 

where the K and F coefficients are defined in table I. The axes and 
the sign conventions employed are shown in figure 1. On applying the 
Laplace transformation 

F(s) = r f(t)e“ s ^dt 
J n 

m 


0 


(where s = a + ico) to both sides of equations (Al), (A2), and (A3) and 
assum in g the In itial conditions to he zero, there is obtained 
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s + Kp 

-k 2 • 

s 


p(s) 


*i8 r (s) 

k 3 

s 2 + K^s 

-V 2 - V 


0(s) 

> = < 

F 2 5 r (s) > 

_-K ? 

-Kqs 2 - K^s 

p 

8 + K 10 s 


t(s) 
1 J 


^(s) 




(A5) 

(a 6) 

(A7) 


(Transfer Functions 

By solving the three simultaneous equations (A5), (A6), and (A7) 
for p(s) there is obtained 


(C 0 's5 + C-^'s^ + Cg's^ + C^'s 2 + C^'sj p(s) = 

(Cj’s 1 *- + C 6 ’s 3 + C 7 's 2 + Cq's^S^s) (A8) 


where the C n ' coefficients are defined in table II . Dividing equa- 
tion (A8) through by Cq' and using the unprimed coefficients C n to 


_ C 1* 


designate the ratio of the primed coefficient to Cq’ (e.g., Cp = 

>\ V 0 

results in the following transfer function for p due to 


and Go = HiL 
e - n t 
°0 


a rudder-deflection input Bp: 


_L = 

&r 


s 


C^s5 + Cgs 2 + CyB + Cq 
^ + Cj_s3 + CgS 2 + C^s + Cj^_ 


(A9) 


Determined in a similar manner, the transfer function for 0 is 



C9S 2 + C 10 s + Cpp 
s^ + Cj_s5 + C 2 s 2 + C^s + C]^ 


and the transfer function for is 


± = 
5r 


Ci pB^ -f- Cp^B 2 + CjjjS 4- Cpr^ 
s ^ + C lS 5 + C 2 S 2 + C-^s + C]^ 


(A10) 


(All) 
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The lateral -acceleration transfer function may "be developed from the 
sideslip transfer function in the following manner: The lateral accelera- 

tion is usually defined as 


ay = v(d? + Dv|r - K^0) 

Substituting equation (Al) into (A12) yields 

ay = W^t) - VK X P 

Therefore, 


or 


fZ = VFp - 
5r 


P(s) 

&r ) 


VKp 


= c l6 s4 


+ + Cjjjs 2 + c 19 s + 

+ Cjs 3 + CgS 2 + C^s + 


c 20 


(A12) 


(A13) 


(A14) 


(A15) 


Modes 

The characteristic equation of lateral motion is obtained by 
expanding the determinant of the coefficients of equations (Al) , (A2) , 
and (A3) and is 


s ^ + 0^5 + C 2 s 2 + C^s + = 0 (Al6) 

This equation can be factored as follows: 

s ^s + 7V]^s + ^2)( s + £ - iT l)( s + i + iq) = 0 (A 17 ) 

where the root s = 0 indicates that the aircraft is insensitive to 
azimuth, the root is the damping -in-roll mode, the root Xg is 

the spiral mode, and the root -g + iq is the Dutch roll mode, an 
oscillatory motion composed of roll, yaw,' and sideslip. The modes of 
this airplane are indicated in table X. 
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An approximate expression for ?v> which gives excellent results is 


_ C4 _ y 


c 3' 


(A18) 


Frequency Response 

The frequency response which is the steady-state response to a 
sinusoidal input consists of an amplitude -ratio and phase-angle relation 
between input and response for various values of cd. It is determined 
by substituting s = ico into equations (A9), (AlO), (All), and (A15), 
respectively, for g, 0, tJc, and ay and is a complex expression. The 

amplitude ratio is the square root of the sum of the squares of the real 
and imaginary parts of the complex expression and the phase angle $ is 
the arc tangent of the ratio of the imaginary part to the real part. 

The amplitude ratio for sideslip angle g is 




( c 8 

8r 

\ 

- 


( c 8 ~ c 6°£ + (p7<p ~ C^ta? 


- C 2 o£ + C]^ 2 + (c 5 o> - Cqo^) 2 

and the phase -angle relation between g and 8^. is 


*08r = tan ' 1 


( V ~ +?p -l / C 3 m - C l°^ \ 

i Cq - Cgo ?J U - C 2 o^ + C J 


The amplitude ratio for roll 0 is 


1 

&r 


fin - c 9^) 2 +( c io^) : 


\ (p - C 2 o£ + C^ 2 +(c^cb - C jcs^Y 

and the phase -angle relation between 0 and 5^ is 

H -tan-! /_£ m \ . tan-1 f By* - By? \ 

?6r l c n - V7 lf» 4 - + o 4 y 


CD' - + C 


(A19) 


(A20) 


(A21) 


(A22) 
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The amplitude ratio for yaw is 


* 


( c 15 “ c l3 a ^) 2 + , c il“ ~ C 12°^) 2 

§r 

\| 

(c lffi 4 - C^a^J 2 + (or 5 - C 2 ou? + C^cd) 2 


and the phase -angle relation between ijr and 5 r is 

t /Ci lim - Ct on>3\ t Ll 


= tan- 1 fa * 0 , - C - tan" 1 ' £ &? . + 

^ ^15 ~ Ci^<d 2 J l C 1 a/ t ' - 


The amplitude ratio for lateral acceleration ay is 


fy 


(° 16 ^ - 

&r 

= \ 



(pl6°^ ~ C l8 a ^ + C 2p) 2 + ( C l9 a> ~ C 17 ^)‘ 


and the phase -angle relation between ay anri 5 r is 


On = tan -1 

X 


( c 19a) - CpyO!? ^ ^ _ 1 j C^O) - CjCD^ 

' j — "fcSJQL j 

^26°^ - C l8 a ^ + C 20/ - C 2 a^ + C 


(A23) 


(A24) 


(A25) 


(A26) 
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appendix b 

DETERMINATION OF LATERAL-STABILITY DERIVATIVES BY USE OF 
FREQUENCY RESPONSES DERIVED FROM CALCULATED TRANSIENT 
RESPONSES TO AILERON DEFLECTION 


Transient responses in roll, yaw, and sideslip to a square-pulse 
aileron deflection were calculated for the airplane whose mass and aero- 
dynamic characteristics are given in table XI. These transients and 
assumed aileron time histories are presented in figure 7. 

The frequency responses D 0/& a , Ity/ S a , a^d P/&a were obtained 

from analysis of the calculated transient responses D f, Bi|r, and p to 
the assumed aileron deflection by use of the equation 



e -stp(t)dt 



e“ s ^q(t)dt 


G(s) 


The Laplace transforms of the responses were obtained by fitting various 
functions, such as polynomials or trigonometric functions, to finite 
sections of the respective output curves and thus performing analytically 
the required integrations from 0 to 00 in finite time intervals as 
described in reference 8. These transforms were then evaluated for a 
number of values of io> (where s = i<r) and frequency-response data 
were obtained. These data have been plotted in figure 8 where they are 
represented by the triangular test points. 

For purposes of comparison, the frequency responses for the airplane 
were calculated from the basic mass and aerodynam i c characteristics and 
are also presented in figure 8 in the form of curves. The agreement 
between these two sets of data is seen to be excellent. 

In addition, the frequency responses were obtained from the calcu- 
lated transients by use of IBM equipment. The Fourier transform of each 
transient was obtained from a numerical evaluation, by means of the 
IBM equipment, of the integral 


e -iojfcF(t)dt 
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between the limits of t = 0 and t = 6.9 seconds to which was added 
an end correction based on an analytical expression of the transients 
from t = 6.9 seconds to infl.nl.ty. A time interval At = 0.05 second 
was used in the numerical integrations. The frequency responses thus 
obtained are shown in figure 8 as square points and appear to be in 
excellent agreement with the frequency responses calculated from the 
basic mass and aerodynamic characteristics and with those derived from 
the transients by the curve -fitting method. 

By use of equations ( 27 ), ( 31 ), (33)* (35)* and (3&) of the body of 
the paper and the frequency-response data derived from the transients by 
both methods, the coefficients Kq_* K^, Kjj., Kg, G 2 * Ky, K^, and K^q 

were obtained, and hence the stability derivatives Gyp* C^p, C^* C j r , 

* *^np* '-'Up* and Cn r - These results are presented in table XII 

along with the values assumed for the calculation of the transient motions. 
Results are also presented for the case where equations (l), ( 2 ), and ( 3 ) 
of the body of the paper were not separated into real and imaginary parts 
but were treated as vector equations and the method of least squares 
applied in a' manner similar to that described in reference 4. The amount 
of work required in the least -squaring process is less when these equa- 
tions Eire treated as vector equations, but a 1 1mi tation is introduced in 
that a parameter such as Kg which is difficult to extract with good 

accuracy cannot be eliminated as it was in the body of the paper when the 
equations were separated. 

A comparison of the values of the derivatives based on frequency 
responses obtained by the curve -fitting method and obtained from the two 
methods of treating the equations of motion indicates very good agreement 
for all the parameters except . The poor agreement for is felt 

not to be significant since C 2 r generally makes a negligible contribu- 
tion to the lateral motions of an airplane and is important only in the 
lightly damped spiral mode of motion which is very difficult to define 
accurately frcan a transient response. 

A si mil ar comparison of the two sets of values based on frequency 
responses obtained with the IBM equipment shows good agreement except for 
the derivatives C 7 , C„ , and Cy . The poor agreement between the 

t r Up 

initially assumed values of C 7 , 

from the frequency responses obtained with the IBM equipment is probably 
due partly to the fact that data were not obtained for frequencies near 
the natural frequency of the airplane with the particular routine used 
for the IBM equipment and partly to smal 1 inaccuracies that exist in the 
data at the frequencies that were included. These small inaccuracies 


'V 


and 




and those calculated 
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undoubtedly could have been reduced somewhat by the use of a smaller time 
interval in the numerical integrations performed by the IBM equipment. 

The use of the frequency-response data which were obtained frcm the 
transient responses by the curve -fitting method and which included fre- 
quencies near the natural frequency of the airplane resulted in good 
agreement for these derivatives. Therefore, it appears that further 
investigation into the effects of the choice of the frequency range to 
be used in extracting stability derivatives frcm frequency-response data 
is warranted. 
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TABLE I.- DEFINITIONS OF THE COEFFICIENTS OF THE 
EQUATIONS OF LATERAL MOTION OF AIRPLANE 


Stability-derivative coefficients : Control-effectiveness coefficients 


Ki 


= % (- 2?) 


K2 = C L h = s/V 


K* = C, 


% 


% =% 
*5 = 

*6 = c i. 


L 2 fc / b ) 2 t2 1 

1 


4-r fc x /b) 2 


K 7 = % 


4t feA) 

t*b 


(V b ) 


2 t 2 


% ~ I xz/ I z 

^ = % Vr (^/ b) 2 


K 10 = c n r 


J|-T 


M”) 


F-, = C Yc . — ■ 
1 1 5 r 2r 


F 2 = C 




2(k v /y 2 t 2 


!Sr p M 

l*b 


F 5 Cn 5 r „ 2 _2 






% = c,6 * 2 yv 


r - r ^b 

^ n 8p « /,. laA 2 „2 


Ja 2 


( k z/ b )' 
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TABLE H.- COEFFICIENTS OF TRANSFER FUNCTIONS 


The equations for the transfer functions are given in appendix A; 
the coefficients C n are obtained from the coefficients C n ' 

by dividing by C 0 ' ^.g., C x = 

C 0 * = 1 - ^Kq 

C l' = K 1 + K4 + K 10 - - KgKg - 

c 2' = k 7 + k 1 k 4- + K 1 K 10 + K k K 10 ~ " K 6 K 9 “ 

C 3 ' = KgK-^ + K]^K»y — + K^Kj^Kjq - K^KgfQ^ — K^K^Ky 

<V = KgKjKiQ - KgKgKj, 

V = Fl (l - K^Kq) 

C 6 ' = F^ + F lKl0 - KqF 2 - F 3 - F-l^ - F lK6 K 8 

V = F lK4 K 10 - F^ + F^ - F 5 K 4 + - F^ 

C 8 ' = F^KgKg + F2K 10 K 2 

c 9 ' = f 2 + f 3 k ^ 

G 10 1 = -^ 2^1 + f i k 5 k 7 ~ f i k 3 + F 3^6 + F 3^i^5 + F 2^lo 
Cu' = F^Nio + Fj^KgKy + F 5 K 5 + F ^ - F^K^ + FjKjKg 

c 12* = F 3 + Ki8 F 2 

Ci 5 « = F 3 K k + F 3 K x - FjK^Kq + F-l^ + FgK^ + F^Kq 

C^' = F 3 K 1 Kj^ - FjKjK^ + FjKjjKj- + F 2 K 1 K^ 

C 15 ' = F 2 K 2 K t + FjKgKj 

c l6* = Co'VFi 

c 17 ' = C 1 *VF 1 - VK l C 5 * 

Ci 8’ = VF-j^Cg ' - VKjCg' 
c 19' = VF^ 1 - VK^C^ 1 

C 20 ' = VF^' - ^03' 
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TABLE HI.- AXRPLANE PARAMETERS USED IN EXAMPLE GIVEN IN BODY OF PAPER 


(a) Known 

M= 0.8 (evaluated for an altitude of 10,000 ft) 
b = 22.6 ft 
S = 130.0 sq ft 
p = O.OOI756 slug/cu ft 
V= 861.74 ft/sec 

I x = 2062 slug -ft 2 


I z = 13,298 slug -ft 2 

m = — = 295*03 slugs 
& 


IXZ = 157 slug-ft 2 

s = ° 


(b) Computed 

T - ik = 1 - 5 se,: 

«> - - 57 - 2 

C L . i = 0.1021 

Kk = = 0.07614 

X X 

KQ = ^ = 0.011806 
I Z 

Kg = | = 0.0374 per sec 


^ - c np 


4-r(k z /b) 


= 0 
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TABLE IV.- FREQUENCY -RESPONSE DATA. FOR EXAMPLE 
GIVEN IN BODY OF PAPER 

(a) Sideslip, {3 


0), 

radians 

sec 


I_ 

Sr 


06/ 

deg 

H - 

§r 

cos 

Bp = 

JL 

§r 

sin 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0.536 
• 574 
.650 

•797 

1.116 

2.105 

4.338 

1.368 

.715 

.464 

-0.1 

-1.1 

- 2.3 

- 4.1 

-8.2 

-20.8 

-109.7 

-158.1 

-166.1 

-169.1 

0.536179 

.574106 

.649929 

.794967 

1.104368 

1.967417 

- 1.462174 

-1.268821 

-.693781 

-.455533 

-O.OOI562 

-.011025 

-.026083 

-.056988 

-.159106 

-.747357 

-4.083763 

-.510100 

-.171676 

-.O8772O 


Cb) Roll, 0 


CD, 

radians 

sec 


1 

Sr 


V 

deg 


1 

Sr 

COB 

B 0 = 

L 

Sr 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

8.761 

4.624 

3.432 

3.111 

3.459 

5-446 

9.719 

2.729 

1.298 

.780 

79-0 

67.7 

57-4 

47.6 

37-1 

18.9 

- 74.6 

- 126.9 

-138.2 

- 143-9 

1.676642 

1.754686 

1.850156 

2.096761 

2.760823 

5.153317 

2.578202 

-1.638908 

-.966715 

-.630341 

8.598903 

4.278325 

2.891147 

2.297575 

2.084203 

1.762687 

-9.370572 

-2.181509 

-.865863 

-.459096 
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TABLE IV.- FREQUENCY -RESPONSE DATA FOR EXAMPLE 
GIVEN IN BODY OF PAPER - Concluded 


(c) Yaw, t 


CD, 

radians 

sec 

1 

JL 

Sr 

1 

*v 

deg 

A|r = 

I_ 

§r 

cos < 5 * 

°r 


JL 

Sr 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

O.223 

.494 

.613 

•773 

1.096 

2.079 

4.299 

1.358 

.710 

.462 

163.4 

171 J -.4 

174.4 
172.8 
168.9 

156.2 

67.3 

18.7 

10.7 
7-1 

-0.213722 

-.491634 

-.610140 

-.767134 

-I.O7525O 

-1.903129 

I.656940 

I.285297 

.698169 

.458498 

0.063714 

.048463 

.059934 

.096775 

.210567 

.837405 

3.966906 

.436504 

.131797 

.057029 


(d) Lateral acceleration, ay- 


03 , 

s. 

Hi/ 

Af , = cos 4 >„ , 

Qrtr 

B„ =/sln V 
'V S3, a y l 

radians 

sec 

Wm 

deg 

ft/sec2 

radian 

ft/sec2 

radian 

1 

107.67 

179.7 

IO7.6736 

0.564215 

2 

121.70 

178.2 

121.6356 

3-928440 

3 

149.83 

176.5 

149.5411 

9.277399 

4 

204.00 

174.0 

202.8943 

21.205800 

5 

322.10 

169.6 

316.7986 

58.143795 

6 

691.42 

156.5 

634.3160 

275.150748 

7 

1628.37 

67.3 

627.9656 

1502.418071 

8 

587.34 

18.6 

556.6625 

187 . 3377 H 

9 

350.66 

10.4 

344.9205 

63.200245 

10 

254.54 

7.2 

252.5535 

31.756598 
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TABLE V.- MODIFIED FORM OP EQUATIONS USED 
HJ EXAMPLE QIVHf HI BOUT OF PAPER 



(-°*P - A + KgBjj) J (Bp) 


} 


^^frAp - K a A <S - “A - A) j 




(27a) 

<26a) 

(29a) 

(30a) 


‘£[mJ 



>- 04 

> 

^«iWj 


s 10 




nF^ - + K 6 a ^ A ji “ K 10 afi fr “ j 



^(«*A ' + V*) j 

^(A) j ( B A + V*)j 

k 3 

^ (V* + Ap*<) , (“AA ' “V 1 ^ j 

J“®L U 






1 

CM 

r^3 

L 

4v|l 

I 2 . 


(31a) 


(32a) 


■a f 

+ “V* - + j 

+ «*V» - (“A ' , 

^ + ^ - V" j (A) j ! 

^ + a?Ajj - B^cA^ + - KjA^ ^ ^ 

|^[H J 

!>' + K^ctA^ + KjlAB^ - <Ab^ ^ ^ 

* 


(3*0 


( 35 a) 


(36a) 
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TABLE VI— LEAST -SQUARES SOLUTION OF EQUATIONS FOR DATA 
OF EXAMPLE GIVEN IN BODY OF PAPER 


(a) Kq, equation (27a) 



-afy + K^Bp = 7.378 



17.562 


I7.562K1 = 7.378 
Kq = 0.420 

o>) *i> equation (28a) 

10 

^^Ap - - oiBp - cnB^) = 1.06l 

0J=1 

lOFq = 1.06l 


F-l = 0.1061 
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TABLE VI.- LEAST -SQUARES SOLUTION OF EQUATIONS FOR DATA 
OF EXAMPLE GIVEN IN BODY OF PAPER - Continued 


(c) % from lateral-acceleration data, equation (29a) 



13,042,324 


2^a^(- TO p) = 5,567,902 


13,042,000% = 5,568,000 
% = 0.427 


(d) F-^ from lateral-acceleration data, equation (30a) 



1.035059 


iof-l = 1.035 


F x = O.IO35 
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TABLE VI.- LEAST-SQUARES SOLUTION OF EQUATIONS FOR DATA 
OF EXAMPLE GIVEN IN BODY OF PAPER - Continued 


(e) Ky and Kq_q, equation ( 31 a) 


21 (- b b ) 2 = 17.562397 

<x=l X PJ 



44.47131 



473.928987 


SZ(“% - %“%)(" Bp)= 856.325838 


o>=l 

10 

2 

10=1 


o% - Kgo%j(o% = 2356.791416 


Kg = 0.011806 = 0 

17.561^ + 44.47K LO = 856.3 
44.471^ + 473.93K 10 = 2356.8 
Ky o 47.443 
K 10 = 0.5217 


(f) F^, equation (32a) 

10F 3 = ^ f Vp + ~ K 10^* " = “ 251 - 


87 


f 5 = -25.187 



4o 
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SUABLE VI.- LEAST-SQUARES SOLUTION OF EQUATIONS FOR DATA 
OF EXAMPLE GIVEN IN BOUT OF PAPER - Continued 


(g) K3, K4, and F 2 , equation (34a) 

10 

X 

( 0=1 


XU 

( B P^r + VV) 2 = 371.762380 


10 

^Z( B f 3 ^ + = -3617.677695 


03=1 

10 
X 

03=1 
LO 


03=1 


+ VV)(" Ajr) = 23.18^790 
51 - oB^) 2 = 35234.438878 

(«VV - ^V)(- A|r) = -209.219323 
(~ Air) 2 = 11.122853 

03=1 V ' 

- k 5 “ 2 [|| 2 )^P b V + Vty) = 33187.734509 
+ °^VV - | ^( cuA ^lr - ^B^v) = -322344.219155 


^5 - K^a?|| A^ = 2422.633922 

371.76238X3 - 3617.67769K4 + 23. 1847952 = 53187.75451 
-3617.67769^ + 35234.43888^ - 209.21932F 2 = -322344.21915 
23.18479K3 - 209 . 21932K^ + 11.12285F 2 = 2422.633922 
K3 = 138.272 K4 = 5.212 


f 2 = 27.636 



IT 
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TABLE VI.- LEAST -SQUARES SOLUTION OF EQUATIONS FOR DATA 
OF EXAMPLE GIVEN IN BODY OF PAPER - Concluded 


(b) Kg, equation (35a) 



8H.5562T 



+ <d 2 A0 - K^ 2 ^ 


+ - K 3 Ap)= 235.117592 


811. 6Kg = 235.1 
Kg = 0.290 


(l) Alternate means of obtaining Kg, equation (36a) 



475.928983 


+ W + K 5^ B V ' 157.774056 


473. 9Kg = 137-77 
Kg = 0.291 
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TABLE VII.- COMPARISON OF COMPUTED AND KNOWN VALUES OF COEFFICIENTS OF 
EQUATIONS OF LATERAL MOTION FOR EXAMPLE GIVEN IN BODY OF PAPER 


Coefficient 

Computed from data 

Known value 

% 

j **0. 42054 
l *0.4269 

| 0.427 

k 3 

138.272145 

138.245 

% 

5.212568 

5.21 

*6 

0.290 

O.3OI7 

k 7 

47.438625 

47.41 

*10 

0.521457 

O.5272 

*1 

10609 
1 *0.1035 

| 0.104 

f 2 

27.636187 

27.65 

f 3 

-25.I87O9I 

-25.22 


a Frcan sideslip data. 

*From lateral-acceleration data. 


TABLE VIII.- COMPARISON OF COMPUTED AND KNOWN STABILITY 
DERIVATIVES FOR EXAMPLE GIVEN IN BODY OF PAPER 


Stability 

derivative 

Computed 

Known 

% 

-1.259 

-1.28 


-0.1490 

-0.149 

° Z P 

-0.4280 

-0.428 

Cz r 

0.02389 

0.0248 

Cnp 

0.3292 

0.329 

°n r 

-0.2732 

-0.279 


O.3183 

0.312 

C »Hr 

0.02979 

0.0298 

Cn 6r 

-0.1748 

-0.175 
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TABLE IX.- COMPARISON OF COMPUTED AND KNOWN COEFFICIENTS OF 
TRANSFER FUNCTIONS FOR EXAMPLE GIVEN IN BODY OF PAPER 


Coefficient 

Computed from data 
( least -squares value) 

Known 

Co’ 

0.999101 

0.999101 

Cl’ 

6.150763 

6.160254 

c 2 ‘ 

50. 934260 

50.972855 

c 3' 

2 53. 483361 

253.2143 

Clt’ 

2 . 193801 

2.190861 

c 5' 

O.IO5995 

O.IO3907 

c 6' 

25.468778 

25.48986 

C 7* 

132.544801 

132.64415 

C8' 

0. 2 67 22 3 

O.2606H 

c 9 ' 

2 5. 718442 

25.729749 

Cio' 

3.636361 

3-952748 

c n' 

-2174.851610 

-2178.769199 

Cl2’ 

-24.86o8l8 

-24.893564 

C 13' 

-136.884812 

-137.264852 

W 

-28.978823 

-30.417543 

C 15* 

-81.654072 

-81.36940 

c l6' 

89.111765 

89.62096 

c lj' 

509.187633 

514.350717 

c l8' 

-4822.026985 

-48l5.4l9l62 

Ci 9 ' 

-26127.620554 

-26138.324577 

c 2 o' 

97.287337 

100.542037 


TABLE X.- COMPUTED LATERAL MODES FOR 
EXAMPLE GIVEN IN BODY OF PAPER 
[Characteristic equation is equation (Al6) of appendix a] 


Damping -in-roll mode, Aq - 5*593 

Spiral mode, -O.OO8668 

Dutch roll mode, -£ ± iq -O.38I + i( 6 . 84 ) 
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TABLE XI.- PARAMETERS OF AIRPLANE USED AS 
EXAMPLE IN APPENDIX B 


Altitude , ft 30,000 

Wing loading, Ib/sq. ft 65 

V, ft/sec 797 

"b, ft 28 

Cl 0.23 

Pb 80.7 

Kx 2 O.OO97 

K Z 2 0.0513 

Kje 0.00145 

Cj^ per radian. - 0.40 

C i T per radian 0.08 

per radian -0.02 

Cn r per radian -0.40 

Cy d per radian -1.0 

P 

C n p per radian 0.25 

C 7 per radian - 0.126 

Cl per radian -0.11 

°a 



TABLE XII.- COMPARISON BETWEEN STABILITY DERIVATIVES USED IN CALCULATION OF 


TRANSIENTS AND THOSE DERIVED FROM TRANSIENTS FOR EXAMPLE IN APPENDIX B 


Derivative 

As Burned value 

Transient analysis , curve -fitting method 

Transient analysis, IEM equipment 

From eqs. of 
present paper 

From eqs. of 
motion as vector eqs. 

From eqs. of 
present paper 

From eqs. of 
motion aB vector eqs. 

c. 


- 0.378 

- 0.397 

-O.388 

- 0.391 

°*r 

0.080 

0.0882 

0.110 

0.130 

0.099 

% 

-0.126 

-O.127 

- 0 . 13 ^ 

-0.129 

-0.128 

°*5 

°a 



-0.100 

-O.098 

- 0.099 


-0.020 

-0.020 

-0.021 

-0.029 

- 0.025 



-0.382 

- 0.397 

-O.786 

- 0.646 

Cup 

0.250 

0.253 

0.257 

0.213 

0.258 

Cy 

Y P 


-1.050 

-1.003 

' -0.600 

1.260 
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Imaginary axis 


D^5 = -a) 2 



5 r = 1 


Eeal axis 


I$B r = -la? 


(a) Vector diagram showing relative magnitudes and positions of input S r 
and output and their derivatives, cd = 0.7 radian/sec. 

Figure 2.- Vector representation of lateral variables. 
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Frequency of rudder-deflection input. 


cd, radiant*/ sec 


Figure 4.- Frequency response in roll due to rudder-deflection input 
for rigid, high-speed airplane of example given in body of paper. 






Frequency of rudder -deflection Input, od, radians /’sec 


Figure 5.- 
rigid, 


Frequency response „ in yaw due to rudder-deflection input for 
high-speed airplane of example given in body of paper. 







Frequency of rudder-deflection, input, m, radians/sec 

Figure 6.- Frequency response in lateral acceleration due to rudder 
deflection input for rigid, high-speed airplane of example given 
in "body of paper. 








Roiling veloci ty y D(f> radians] 





Time, isec 


(c) Angle of sideslip 


Figure 7.- Concluded. 



Amplitude ratio 


rai i to 


Calculated from mass and aerodynamic 
characteristics of airplane 

A Tr ans ient analysis , curve-fitting 

□ Transient analysis, IBM egulpnent 

(At = 0.05 seo) 


cd, radlans/sec 


(a) Amplitude ratio IJ-I. 


Figure 8.- Frequency responses due to aileron-deflection input for the 

airplane described in table XI. 











Calculated from mass and aerodynamic 
characteristics of airplane 

Transient analysis,, curve-fitting 
method 

Transient analysis, IBM equipment 
(At = O.Op Bee) 
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